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Shen-ge  Wang  and  Nicholas  George 
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ABSTRACT 

Automatic  recognition  of  faces  is  now  a  reality,  in  our  experiments  we  seek  to 
establish  bounds  on  the  information  content  required  as  a  function  of  the 
configuration  of  the  "machine."  In  general  with  the  direct  image,  excellent 
recognition  is  obtained  for  information  capacity  below  Harmon's  (1973) 
benchmark  value  of  500.  Our  main  data  base  is  gathered  from  25  people 
consisting  of  three  views  (front,  30°  turn,  60°  turn)  and  10  separate  exposures 
for  each  view.  Typically,  these  are  grouped  into  a  learning  set  of  5  and  a  testing 
set  of  5.  In  the  first  series  of  experiments,  an  optical  Fourier  transform 
configuration  is  used  which  includes  a  ring-wedge  photodetector  interface 
together  with  neural-network  software.  Good  recognition  is  obtained,  and  we 
report  on  the  training  using  front  and  side  views.  In  the  second  series  of 
experiments,  frame-grabbed  images  are  used,  again  with  the  neural  network 
software.  Tests  are  made  of  the  information  capacity  required  as  a  function  of 
the  image  coding.  Comparative  error  rates  are  presented  for  recognition  in  the 
following  cases:  with  256  grey  levels  or  using  global  binarization  or  using  an 
error-diffusion  halftoning  algorithm.  Particularly  good  results  are  obtained  from 
the  error-diffusion  form  of  the  imagery. 
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OPTICAL  FOURIER  TRANSFORMS 
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PARTICULATE  SORTING  USING  MATCHED  FILTERS 


PARTICULATE  SORTING  USING  MATCHED  FILTERS 
Scott  D.  Coston  and  Nicholas  George 

Institute  of  Optics 
University  of  Rochester 
Rochester,  NY  14627 

ABSTRACT 

We  present  theoretical  and  experimental  results  using  matched  filters  to  sort 
between  two  classes  of  spherical  particles.  In  particular,  matched  filters  are 
used  to  distinguish  between  dielectric  and  metal  spheres.  Three  separate  light 
scattering  regions  are  compared  in  order  to  identify  the  best  scattering  angles. 
These  include  the  forward,  the  near  forward,  and  the  backscatter  regions. 
Scattering  theory  is  used  to  calculate  the  correlation  and  cross-correlation 
peaks  for  the  pertinent  cases  in  the  output  plane  of  the  matched  filter.  An 
interesting  finding  that  has  not  appeared  in  the  literature  is  that  the  forward 
scattering  region  is  an  ineffective  area  in  which  to  use  a  matched  filter  for 
distinguishing  between  these  particles.  Of  the  three  scattering  regions,  the 
backscatter  is  shown  to  be  the  most  effective.  Experimental  confirmations  of 
these  findings  have  been  made  using  holographic  recorded  matched  filters  in  a 
standard  Fourier  transform  configuration.  These  will  be  described,  including  the 
successful  matched  filter  discrimination  utilizing  the  backscattered  radiation. 


PARTICULATE  SORTING  USING 


Presnted  at  the  1991  Optical  Society  of  America  conference  In  San  Jose. 
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SCATTERING  GEOMETRY 
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MIE  SCATTERING  FORMULAS 


The  scattered  fields  in  the  far  zone  can  be  calculated  using  the  Mie  formulas1.  The  coefficients  an  and  bn  contain  only  information  about  the  sphere, 
while  the  functions  Tn  and  ixn  contain  the  angular  information.  Calculation  of  the  coefficients  are  accomplished  using  recursion  formulas2  for  an,  bn, 
V  and  «n. 

1 . )  G.  Mie,  "Beitrage  zur  Optik  Truber  Medien,  Speziell  Koloidaler  Metalldsungen,"  Ann.  Physik  25, 377-445,  (1908). 

2. )  W.  D.  Ross,  "Computation  of  Bessel  Functions  in  Light  Scattering  Studies,"  Appl.  Opt.  11, 1919  (1972). 


ELECTRIC  FIELD  MAGNITUDE 


This  figure  shows  the  log  of  the  magnitude  of  the  x-polarized  scattered  field.  The  solid  curves  was  calculated  for  a  dielectric  sphere  with  index  of 
refraction  1.5.  The  dashed  curve  is  for  a  perfectly  conducting  sphere.  Both  spheres  have  a  radius  of  10  |im.  Note  the  similarity  of  the  two  plots  near 
the  forward  scattering  region.  Past  6=5°,  however,  the  scattering  pattern  for  the  metal  sphere  decays  to  a  constant  value,  proportional  to  the  radius, 
while  the  scattering  pattern  for  the  dielectric  sphere  continues  to  oscillate  throughout  the  region  shown. 


ELECTRIC  FIELD  PHASE 


ELECTRIC  FIELD  PHASE 


In  this  figure,  we  show  the  phase  of  the  scattered  field  polarized  in  the  x-direction  for  the  dielectric  sphere.  The  phase  begins  at  approximately  -1.7, 
decreases  slightly,  and  then  oscillates  between  -n  and  +7t.  Note  also  the  relatively  constant  period  of  the  phase. 


MAKING  THE  MATCHED  FILTER 


■I 


To  include  both  the  magnitude  and  phase  of  the  backscattered  scattered  field  in  a  sorting  routine,  we  propose  to  use  a  matched  filter.  The  construction  of 
a  matched  filter  is  shown  in  this  figure.  The  backscattered  field,  at  the  film  location,  is  denoted  by  Esl.  A  reference  beam,  propagating  at  an  angle  to 
the  incident  beam,  is  interfered  with  Egj.  The  transform  lens  is  used  to  produce,  at  the  back  focal  plane,  the  optical  transform  pattern  of  the  reference 


INTENSITY  VALUE  AT  THE  CORRELATION  PLANE 


Mathematically,  the  density,  D,  of  the  film  after  developing,  will  be  given  by  the  intensity  of  the  two  interfering  beams,  I.  The  third  term  in  the 
expression  for  I  constitutes  the  matched  filter.  If  the  reference  beam  is  removed,  and  a  second  sphere  with  scattered  field  Eg2  is  examined,  the  output  at 
the  back  side  of  the  film  will  be  given  by  Es2D.  The  third  term  in  the  quantity  Es2D  is  a  reconstructed  reference  beam,  modulated  by  the  cross 
multiplication  of  the  two  scattered  field,  Esl  and  Es2.  The  transform  lens  will  therefore  produce  the  optical  transform  intensity  of  this  quantity.  The 

intensity  of  the  optical  transform  pattern,  at  the  correlation  plane,  is  Isl.s2-  The  symbol  ®  represents  a  correlation  and  f  represent  a  Fourier  transform. 


center  with  no  surrounding  structure. 


BACKSCATTER 
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whereas  jagged  spikes  represent  a  dielectric  sphere. 


BACKSCATTER 


This  plot  shows  Idiel-DIEL  usinS  the  backscattered  field  of  a  dielectric  sphere  versus  the  normalized  transform  frequency,  fc.  As  with  the  metal 
sphere  previously,  there  is  a  sharp  peak  at  the  center.  Therefore,  the  dielectric  sphere  can  also  be  used  to  make  the  matched  filter.  Analysis  of  the 
intensity  at  the  correlation  plane  is  the  :>ame  as  for  the  matched  filter  made  with  the  metal  sphere. 
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RECOVERY  OF  PARTICLE  SIZE  DISTRIBUTIONS  BY 
INVERSION  OF  THE  OPTICAL  TRANSFORM  PATTERN 
Scott  D.  Coston  and  Nicholas  George 

Institute  of  Optics 
University  of  Rochester 
Rochester,  NY  14627 

ABSTRACT 

Experiments  are  presented  for  recovering  the  histograms  of  particle  size  from 
the  optical  transform  intensity  of  spherical  particles.  We  describe  a  novel 
inversion  formula,  and  we  provide  computational  details  for  obtaining 
particulate  sizes  from  the  transform  intensity.  The  intensity  is  given  by  an 
integral  over  the  size  distribution  multiplied  by  an  Airy  function  (transform 
intensity  for  a  single  particle).  The  inversion  formula  is  derived  from  the 
intensity  using  Bessel-delta  identities  and  Abel  transforms.  By  multiplying  the 
intensity  by  an  appropriate  kernel  and  integrating  over  the  diffraction  angle, 
recovery  of  the  particle  size  distribution  is  obtained.  In  the  experiments,  we  use 
chrome  masks  with  large  numbers  of  tiny  circular  apertures  or  discs.  These 
have  been  generated  with  three  separate  distributions:  the  gamma,  the  log¬ 
normal,  and  the  bi-modal  exponential.  A  discussion  will  be  made  of  the 
calibration  procedure  for  the  ring-wedge  detector  as  well  as  limitations  due  to 
particle  size  and  number.  Results  using  the  inversion  formula  show  excellent 
agreement  with  the  known  values  of  the  size  distributions. 


RECOVERY  OF  PARTICLE  SIZE  DISTRIBUTION 


Presnted  at  the  1991  Optical  Society  of  America  conference  In  San  Jose. 


INCIDENT  PARTICLES  TRANSFORM  RWD 

BEAM  P(a)  LENS  1(0) 


This  figure  shows  the  general  geometry  for  recording  the  optical  transform  pattern  for  a  collection  of  particles,  incident  on  by  a  plane  wave.  The 
particles  are  described  by  the  size  distribution  p(a).  The  detector  records  the  optical  transform  pattern  of  the  particles,  produced  at  the  back  focal  plane 
of  the  transform  lens.  The  scattering  angle  8  is  measured  from  the  optical  axis  to  a  point  on  the  detector.  The  detector  consist  of  32  semi-angular  rings 
and  32  wedges.  However,  since  the  transform  pattern  for  a  spherical  particle  is  circularly  symmetric,  only  the  ring  detectors  are  used. 


DERIVATION  OF  THE  INVERSION  RESULT 


O  ^ 

T* 

«  o> 
c  o> 
o>C- 

to 

E  g 


O  Q. 

•  o 
S3  . 

-  Q. 

c  a 

i: 

o 

CD  w 

a  S 


?! 


2:  <r> 
jc  o 
(0  o 


illto 

§  1 1  i  I 

2  £"<£  .c  > 

«  .S  «  o  C 

c  S  ^  c  S 

£5  ca  a  O 
c  <+M  c  o  o 

O  ®  |  9  2 

£  £*  ^  2  c 

11  &i 

1  1 1 1 S 

fc*  O  *> 

■g  -S  s 

5  g>“  «- 
•“  O  GO  S 

O  *0  7?  c 
.  o  .2  *r  O 

m  L*  G 

C  4)  O 

o  .9*  5  S  -2 
>  5  &  «  3 

3  g  J  "  o 

«  C  ±1  O  “5 

a.  »<*-•-  u 

S  os  o  .S 

-•  &&.S2 
e  e  £  c  S 

'M  *8  8  ”  0 
•a  »  u  £>  5 
c  .o  «-.  >  <2 

8  'fi  s  ®  -S 

**  =3  ~  8  © 

c  g  ~  3  £ 

co  .H  is  ^3  .2 

~  ««  3  >  . 

2  a>  «o  d 

(Tj  q  g  crt  4) 

1  If  2  |i> 

J  op  °  « 

o  «  >,  —  S  St; 

2  •S  .t:  .£  £  o 

O  h  C  O  u  3 

•2  «  §  £  5  o 
C  O  S  -  S| 

-  6)  !  1  ®  O 

v  o'  o  ’S1J2  c 

V5  £5  2P  fan  S 

ft  C  oo  3  c  -2 

<u  o  ts  ^ 

QQ  C  4>  .C  **  O 

o  §  ■£  *  %  § 

H  cc  _,  %  O  Vt 
*~  -S  ~  CO 

>r  o  i  ®  s  § 

els  MS  O 

S  /S  M  c  _  -c 
-«©  e  *«  3 

ig  >*.S  -S  «  s 
;c  -°  «  5  c  ” 

’♦*  *0  *52  m  o 

c  S  *5  «j  d  u 
£  §  q  -c  c 

a"  g  c  o  bo-3 
,jC  o  .2  b.  '*g  § 

H  g.<2 s  ll 

g.  S'  o  g  c 

P€lhs 

c  c  <s  o  w*  52 

o  «  ®  c  ^  O 

w  °  ^  ‘c  ^ 

‘  cO  M  /-v  e  ,5  t/) 

m  O  A  3  w*  -s 

2Sv®£5 
g  S  5  ^  o  S 
S  -o  i*  o  S 

«2  o  -T  '55  'C 

liS  I®  i 

fi  .o  h-P.£  «  .5 


INVERSION  -  MULTIPLE  DISCRETE  SIZES 
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INVERSION  RESOLUTION 
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involved  in  calculating  the  inversion  integral. 


RECOVERY  OF  GROSS  STRUCTURE  -  SMALL  0 


To  verify  the  results  in  the  resolution  plot,  we  created  a  particle  size  distribution  consisting  of  500  discrete  sizes.  Using  8Max=0-05  rads  (5.0  pm 
resolution),  the  resulting  recovered  distribution  was  calculated  as  shown  in  this  plot.  The  ensemble  average  of  this  500  discrete  distribution  is  also 
plotted.  Note  the  smooth  nature  of  the  recovered  distributions.  The  following  plots  show  the  effects  of  increasing  the  value  of  0MAX. 


DISTRIBUTION  RECOVERY 


As  before,  but  with  Bmax1*0*2  rads‘  1116  rcsolution  is  now  10  ^m-  Note  the  beginnings  of  some  extra  detail  in  the  recovery. 


RECOVERY  OF  FINE  STRUCTURE  -  LARGE  6 


(e)d 


As  before,  but  with  0jyfAX=l.O  rads.  The  resolution  is  now  0.2  |im.  Note  the  jagged  appearance  of  the  recovered  distribution, 
changes  in  the  distribution  represent  the  discrete  nature  of  the  particle  size  distribution.  Examination  of  the  actual  particle  sizes  us< 
distribution  reveafcthat  this  fine  scale  structure  should  be  present. 


KAKIIULt  MAS* 


aperture  on  a  chrome  coated  glass  slide  and  tiny  chrome  dots  on  a  gl 
example.  The  area  containing  the  particles  measures  2.5  mm  square. 


INVERSION  -  OPTICAL  DATA  INPUT 


In  the  first  experiment,  we  used  a  log-normal  distribution  consisting  of  6000  particles.  Using  0MAX=O.O5  radians,  the  recovered  distribution  versus 
particle  size  is  shown  in  this  plot.  In  addition,  we  show  the  ensemble  average  curve  for  the  log-normal  distribution.  Note  the  excellent  agreement 
between  the  two  curves. 


INVERSION  -  OPTICAL  DATA  INPUT 
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In  a  second  experiment,  we  varied  the  number  density  of  a  log-normal  distribution  to  confirm  the  linearity  of  the  inversion  formula.  In  this  figure  we 
plotted  the  recovered  distributions  versus  particle  size  using  Omax^-S  radians.  Note  the  similarity  of  the  four  plot,  differing  only  in  the  height. 


INVERSION  RESULTS 


Here  we  show  a  comparison  of  the  known  particle  numbers  used  in  the  previous  figure,  to  number  obtained  by  integrating  the  recovered  distribution 
over  all  particle  sizes.  We  see  that  there  is  excellent  recovery  of  the  particle  number. 


EVALUATION  OF  PHOTOLITHOGRAPHY 


In  a  final  experiment,  we  examined  the  ability  to  distinguish  small  changes  in  the  particle  size  distribution.  For  this  experiment,  we  made  masks  with 
apertures  and  chrome  dots  using  the  same  original  masters.  However,  due  to  variation  in  the  construction  process,  slight  differences  in  the  sizes  of  the 
particles  can  be  measured.  Under  microscopic  examination,  the  discs  were  consistently  5%  smaller  in  size  than  the  corresponding  apertures.  In  this  plot 
we  show  the  recovered  distribution  using  the  aperture  and  disc  mask.  The  value  of  0MAx=0  045  radians  was  used.  The  discs  are  approximately  5% 
smaller  in  radii  than  the  apertures. 


BEAM  SPLITTER  CUBE  FOR  WHITE-LIGHT  INTERFEROMETRY 
Keith  B.  Farr  and  Nicholas  George 

Institute  of  Optics 
University  of  Rochester 
Rochester,  NY  14627 

ABSTRACT 

We  describe  a  technique  for  constructing  beam  splitter  cubes  in  which  we 
monitor  the  optical  path  difference  with  a  wavelength  tunable  laser.  A  white- 
light  interferometer  should  have  path  lengths  equal  to  within  X/4.  A  larger  path 
difference  results  in  reduced  fringe  visibility  and  loss  of  signal-to-noise  ratio. 
The  beam  splitter  must  be  manufactured  such  that  the  two  paths  through  glass 
are  equal.  We  discuss  the  tolerance  placed  on  the  cube.  One  can  assemble 
the  cube  from  two  right  angle  prisms  in  optical  contact  and  in  near-perfect 
alignment.  Alternatively,  we  adjust  the  OPD  to  zero  by  sliding  the  prisms  into 
proper  alignment  on  a  layer  of  uncured  cement.  Using  a  dye  laser,  we  have 
readily  achieved  acceptable  path  differences.  This  method  of  measuring  OPD 
with  a  tunable  laser  is  shown  to  be  useful  for  alignment  of  interferometers  and 
measurement  of  aberrations  in  lens  systems.  A  laser  diode  can  be  tuned  by 
changing  its  operating  current  or  temperature.  We  report  results  of  optical  path 
difference  measurements  using  a  laser  diode. 
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Scope 


Dye  laser  technique  for  measurement  of 
interferometer  path  difference 


Calculation  of  path  difference 


Beam  splitter  geometry 


Extra  geometrical  length  in  path  II: 
h=OB-OA=CD 


By  examination  of  triangle  O'CD  we  find 


_  A+e 

~^r 


Beam  splitter  assembly 


Fringes  observed  on  a  reticule 


Visually  measure: 


h  <  0.25  \xm 


Experimental  Results 

Beam  splitter  path  difference  measurements 


everse 


wmie-iigni  imeneromeiry: 

Spatial  transform  system 


If  interferometer  paths  are  equal 


Summary 


Fig.  1  The  path  difference  h  in  a  two-beam  interferometer  can  be  measured  using  a 
wavelength-tunable  laser.  A  tilt  6  is  introduced  to  cause  a  spatial  fringe  pattern  to  be 
formed. 

Fig.  2  When  the  wavelength  is  increased  the  phase  decreases  monotonically  causing 
the  fringes  to  move.  With  each  2jc  phase  change  the  pattern  moves  one  fringe  width. 
The  fringe  count  is  proportional  to  the  geometrical  path  difference  h. 

Fig.  3  A  beam  splitter  cube  is  constructed  of  two  identical  right  angle  prisms.  If  they 
are  placed  into  contact  to  form  a  perfect  cube  then  the  path  difference  will  be  zero. 

Two  possible  misalignment  are  a  cement  layer  of  thickness  A  and  slide  misalignment  e 
defined  by  the  distance  between  the  altitudes  of  the  two  prisms.  Even  if  a  cement  layer 
is  used  one  can  always  find  a  position  at  which  the  path  difference  is  zero. 

Fig.  4  The  beam  splitter  cube  is  assembled  under  interferometric  control.  Two  prisms 
are  assembled  with  a  layer  of  uncured  cement  and  placed  into  a  special  jig  which  fixes 
one  prism  while  allowing  the  other  to  be  moved  by  a  differential  micrometer.  A  dye 
laser  illuminates  the  cube  and  an  interference  pattern  is  formed  by  the  light  reflected 
from  the  cube  faces.  As  the  laser  wavelength  increases  the  fringe  pattern  moves  and 
the  fringes  are  counted.  The  prisms  are  then  slid  into  a  position  where  the  fringes  are 
stationary.  At  this  point  the  path  difference  is  zero.  The  cement  is  set  using  an 
ultraviolet  lamp. 

Fig.  5  Experimental  results  of  path  difference  measurements  of  several  beam  splitters. 

Fig  6.  The  spatial  transform  system  is  an  example  of  a  white-light  interferometer.  The 
pinhole  located  a  distance  h  from  the  optic  axis  is  illuminated  by  a  source  of  bandwidth 
AX  centered  at  Xo-  The  interferometer  consists  of  a  beamsplitter  and  two  right-angle 

prisms  tuned  at  90^  with  respect  to  each  other.  With  this  arrangement  the  point  source 
at  h  appears  as  two  point  sources  separated  by  2h  and  which  interfere  coherently. 

The  fringe  achromat  is  designed  such  that  the  fringes  at  the  detector  do  not  change 
with  wavelength.  For  this  reason  the  spatial  transform  system  is  particularly  good  for 
measurements  of  visibility.  Such  an  interferometer  requires  that  the  path  difference  in 
the  beamsplitter  be  from  1  to  10  pm. 


BACKSCATTER  FROM  A  TILTED  ROUGH  DISC 


BACKSCATTER  FROM  A  TILTED  ROUGH  DISC 
Donald  J.  Schertler  and  Nicholas  George 

Institute  of  Optics 
University  of  Rochester 
Rochester,  NY  14627 

ABSTRACT 

We  present  theoretical  and  experimental  results  for  the  backscattered  intensity 
from  a  roughened  metallic  disc  tilted  at  an  arbitrary  angle  to  the  incident 
illumination.  As  the  rms  roughness  is  increased  from  1pm  to  10pm  the  zero-tilt 
intensity  decreases  by  more  than  1000  but  at  a  75°  tilt  it  is  shown  to  increase  by 
more  than  100.  We  also  investigate  the  decorrelation  of  speckle  with 
wavelength  at  various  tilt  angles.  Tilting  the  surface  to  the  illumination  and  then 
observing  the  intensity  pattern  in  the  specular  direction  gives  a  means  of 
characterizing  the  roughness  parameters.  In  performing  the  calculation  the 
effect  of  shadowing  must  be  included.  We  analyze  this  effect  by  comparing  the 
calculated  results  with  shadowing  to  those  without  shadowing.  The  effect 
becomes  significant  for  moderate  tilt  angles  (>60°)  and  large  roughness.  In 
addition  we  review  recent  developments  in  enhanced  backscatter  as  applied  to 
rough  surfaces.  Several  techniques  for  producing  the  roughened  surfaces  are 
also  described  including:  (1)  beadblasting,  (2)  photolithography/ion-etching  and 
(3)  exposing  of  photo-resist  to  controlled  speckle  patterns.  These  techniques 
give  us  the  ability  to  produce  surfaces  with  rms  roughnesses  ranging  from  1pm 
to  10pm  and  correlation  lengths  from  1pm  to  50pm. 
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The  Scattered  Electric  Field  Component 


The  Slope-Height  Term: 


Gaussian:  r(p)  =  exp(-p2/lc) 

Exponential:  r(p)  =  exp(-|pj/lc) 


Result  for  aBS  at  0,=O 

Gaussian  surface: 

7ia2 

-fcj 

Exponential  surface: 

°BS 

=2f2adeYf x ) 

Tia2 

l  X  )  V 2no) 

I2 

=  2  c 

kV 

The  cross  section  of  the  Gaussian  surface  at  normal 
incidence  is  much  larger  than  that  of  the  exponential 
surface. 

The  cross  section  of  the  exponential  surface  falls  off 
much  more  rapidly  with  increasing  roughness,  o. 

For  non-normal  incidence  we  must  consider  the  effects 
of  shadowing. 


Shadowing  of  a  surface  is  a  function  of  the  product  tan0, 
times  the  rms  slope  of  the  surfacet 


The  rms  slope  is  given  by 


For  the  Gaussian :  =  V2  <j/lc 


For  the  exponential  surface  this  formula  gives 


/  dVl 
l  dp2J 


S(p)  -  j-  exp( — |p|  /  lc ) 


The  rms  slope  can  be  determined  from  computer 
simulation  and  approximation  of  the  delta  function.  We 
have  found: 

The  rms  slope  for  the  Exponential :  =2.8  a/ 

tP.  Beckmann,  IEEE  Trans.  Antennas  &  Propag.  AP-13  (1965). 

B.  G.  Smith,  IEEE  Trans.  Antennas  &  Propag.  AP-15  (1967). 


The  Shadowing  function  is  the  fraction  of  the  surface 
that  is  illuminated  at  a  given  incident  angle. 


Backscattering  Cross  Section  versus  Tilt  Angle 
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Measured  Backscattering  Cross  Section  vs  Tilt  Angle 
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Page  2.  Rotating  Discs 

This  series  of  pictures  illustrates  the  problem.  Shown  are  three  discs  tilted  at  various 
angles  tothe  illumination.  The  discs  are  viewed  through  a  beam  splitter  so  that  we  see  the  direct 
backscatter.  The  tilt  angles  are  a)  zero  degrees,  b)  20°,  c)  40°  and  d)  70°.  The  large  disc  has 
a  diffuse  white  surface.  The  smaller  two  are  one  inch  diameter  silver  discs  that  have  been 
roughened:  the  upper  left  at  2pm  rms  roughness  using  a  bead-blasting  process,  and  the  lower 
at  8pm  using  ion-beam  etching  of  a  speckle  pattern  exposed  in  photo-resist.  As  the  discs  are 
tilted  the  2pm  disc  scatters  less  and  less  toward  the  illuminating  direction  while  the  8pm  and 
diffuse  white  discs  continue  scattering  strongly.  At  70°  the  8pm  disc  is  now  considerably 
brighter  than  the  othertwo.  This  is  explained  in  terms  of  the  roughness  parameters  of  the  discs. 


Page  3.  Vector  Potential 

In  our  analysis  the  disc  is  centered  at  the  origin  of  our  coordinate  system  in  the  xy  plane. 
It  is  illuminated  with  a  y-polarized  plane  wave  whose  k  vector  lies  in  the  xz  plane  and  forms  an 
angle  Q.  with  the  z-axis.  The  backscattering  cross  section  is  defined  as  shown  where  IBS  is  the 
backscattered  intensity  at  P  a  distance  R0  away. 

We  use  a  vector  potential  analysis  so  that  we  can  write  H  as  the  curl  of  the  vector 
potential.  At  P  the  vector  potential  is  an  integral  of  the  surface  current  density,  K.  We  use  these 
to  derive  an  expression  for  the  far  zone  electric  field. 

Page  4.  Far  Zone  Electric  Field  Component 

The  far  zone  y-component  of  the  scattered  electric  field  is  given.  The  integral  is  over  the 
illuminated  surface  of  the  disc  and  contains  terms  involving  the  slope  and  height  of  the  surface 
at  the  point  of  integration. 

The  intensity  is  found  by  taking  an  ensemble  average  of  the  magnitude  squared  of  the 
field.  We  have  grouped  the  terms  involving  the  surface  height  and  slope  into  the  function,  F. 


Page  5.  The  Slope-Height  Term 

The  function  F  is  the  ensemble  average  of  the  slope  and  height  terms  at  two  surface 
points.  It  can  be  expanded  to  three  terms  as  shown.  The  cos2  factor  is  the  result  for  a  smooth 
surface,  i.e.,  h=0.  The  average  of  the  exponential  term  is  the  joint  characteristic  function  of  the 
two  heights,  h’  and  h".  The  other  terms  affect  the  result  for  non-normal  incidence  due  to  the 
sine  term.  The  first  of  these  two  terms  contains  averages  of  the  slope  at  a  given  point  times 
the  phase  difference  between  two  surface  points.  The  last  term  correlates  the  slope  at  two 
points  with  this  phase  difference. 

To  fully  evaluate  this  expression  we  require  a  surface  correlation  function. 


Page  6.  Surface  Correlation 

The  surface  correlation  function  is  given  by  the  normalized  ensemble  average  of  the 
heights  where  a  is  the  rms  height  variation  and  p  is  the  distance  between  h'  and  h".  We  have 
assumed  that  h  is  characterized  by  a  stationary  Gaussian  random  process.  The  correlation 
length  is  denoted  by  lc.  We  have  examined  two  types  of  surface  correlations:  a  Gaussian  and 
an  exponential.  The  graph  illustrates  the  difference  between  these  surfaces.  The  Gaussian 
has  slowly  varying  and  continuous  slopes  while  the  exponential  has  steep  and  discontinuous 
slopes. 


Page  7.  On-axis  Results 

The  on-axis  results  for  the  backscattering  cross  section  for  the  two  surfaces  are  given. 
These  hold  for  o»\.  We  note  that  the  result  for  the  Gaussian  surface  is  independent  of 
wavelength  while  the  exponential  result  varies  as  X2.  Also,  the  Gaussian  result  is  much  larger 
than  that  of  the  exponential.  The  cross  section  for  the  exponential  surface  falls  off  rapidly  with 
increasing  roughness,  o. 

To  determine  the  off-axis  cross  section  we  must  consider  the  effect  of  shadowing. 


Page  8.  Shadowing 

The  shadowing  of  a  surface  is  a  function  of  the  product  tan9j  with  the  rms  slope  of  the 
surface.  The  rms  slope  is  proportional  to  the  square  root  of  the  second  derivative  of  the 
correlation  function.  For  the  Gaussian  surface  this  evaluates  to  the  result  shown.  For  the 
exponential  surface,  however,  the  second  derivative  of  the  correlation  function  gives  a  delta 
function  at  the  origin.  The  rms  slope  can  be  evaluated  by  simulating  exponential  surfaces  on 
a  computer  and  approximating  the  delta  function  based  on  minimum  resolvable  surface 
features.  We  have  found  that  the  rms  slope  forthe  exponential  surface  is  given  approximately 
by  the  result  shown.  We  see  that  this  value  is  much  larger  than  that  for  the  Gaussian  surface 
as  was  demonstrated  on  page  6. 


Page  9.  Shadowing  continued 

These  results  are  used  to  determine  a  shadowing  function  which  is  the  fraction  of  the 
surface  which  is  illuminated  at  a  given  incident  angle.  Shown  in  the  top  graph  is  a  simulated 
surface  illuminated  at  0r  Below  this  is  a  function  indicating  by  the  value  of  1  those  portions  which 
are  illuminated  and  by  0  those  which  are  in  shadow.  Thus  ~50%  of  the  surface  is  illuminated 
at  this  angle.  Carrying  this  out  for  other  angles  and  over  many  surfaces  we  obtain  an  average 
shadowing  function  S(0)  shown  in  the  lower  graph  forthe  rms  roughness  and  correlation  length 
given.  The  exponential  surface  is  represented  by  circles  and  the  Gaussian  by  squares.  For 
the  same  roughness  parameters  significant  shadowing  occurs  at  a  much  smaller  angle  forthe 
exponential  surface  than  forthe  Gaussian.  We  use  this  shadowing  function  in  the  result  forthe 
backscattering  cross  section  for  large  incident  angles. 


Page  10.  The  Backscattering  Cross  Section 

The  results  for  the  backscattering  cross  section  versus  tilt  or  incident  angle  are  shown 
for  both  types  of  surfaces.  Each  is  given  for  a  fixed  correlation  length  and  for  four  values  of  the 
rms  height.  The  normalized  cross  section  is  plotted  on  a  logarithmic  scale.  The  Gaussian  cross 
sections  fall  off  with  increasing  tilt  angle,  yet  requiring  a  larger  tilt  angle  for  the  fall-off  to  occur 
with  increasing  roughness.  For  the  ranges  shown  no  significant  shadowing  occurs. 

For  the  exponential  surface  we  have  included  two  sets  of  results:  those  where  we  have 
included  shadowing  in  the  calculation,  signified  by  the  solid  lines  and  those  without  shadowing 
indicated  by  the  dashed  lines.  The  cross  section  remains  significant  for  all  the  roughnesses 
shown  at  nearly  all  angles  of  illumination.  For  large  roughnesses  the  cross  section  actually 
increases  with  tilt  angle. 


Page  1 1 .  Measured  Results 

Shown  here  are  the  measured  backscattering  cross  sections  of  the  three  discs  shown 
on  Page  2.  The  diffuse  white  surface  exhibits  the  cos2  behavior  of  a  blackbody  scatterer.  The 
2pm  disc  shows  a  more  rapid  fall  off  yet  still  has  a  substantial  value  at  large  angles.  The  8pm 
disc  follows  the  diffuse  white  disc  up  to  about  60°  where  it  begins  to  increase  up  to  ~70°  and 
then  falls  off  again.  The  large  cross  section  at  0°  and  the  increase  at  large  angles  indicates  that 
this  surface  may  be  a  composite  surface  containing  both  Gaussian  and  exponential  correlation 
components. 
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ABSTRACT 

We  present  analysis  for  the  recovery  of  intensity  images  degraded  by  an 
incoherent  point  spread  function  (psf)  consisting  of  a  set  of  randomly  distributed 
impulses.  In  initial  experiments  with  as  many  as  100  irregularly  spaced 
impulses  we  obtain  excellent  recovery  of  the  original  object  provided  we  have 
precise  knowledge  of  the  psf.  In  practice,  there  is  little  or  no  difficulty  presented 
by  the  zeros  of  the  multiple  point  blur  transfer  function.  We  describe  some 
interesting  methods  of  counting  these  zeros  over  an  appropriate  region  of  the 
transform  plane.  In  the  second  portion  of  our  presentation  particular  emphasis  is 
placed  on  the  effects  of  errors  made  in  the  determination  of  the  point  spread 
function.  We  simulate  measurement  errors  in  the  recording  of  the  psf  by 
randomly  perturbing  the  location  of  each  impulse  about  the  actual  position  used 
to  create  the  blurred  image.  Analysis  and  reconstructions  based  on  variations  in 
the  measured  psf  are  presented.  Examples  include  the  statistical  deviation  of 
the  measured  impulses  from  their  true  locations  and  the  number  of  impulses 
present  in  the  degradation. 


Image  Deblurring  for 
Multiple-Point  Impulse 
Responses* 


Bryan  J.  Stossel 
Nicholas  George 


Institute  of  Optics 
University  of  Rochester 


^This  Research  is  Supported  by  the  US  Army 

Research  Office 


Presented  at  the  1991  annual  meeting  of  the  Optical  Society  of  America 

San  Jose,  Ca. 


Imaging  Model: 


OBJECT  SPACE: 

g(x,y)  =  f(x,y)*h(x,y)  +  n(x,y) 
where 

g(x,y)  =  observed  image 
f(x,y)  =  original  object  scene 
h(x,y)  =  impulse  response 
n(x,y)  =  noise 

FREQUENCY  SPACE: 

G(Vv)  =  F(My)  H(fx>fy)  +  N(fx,fy) 
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Charcterization  of  the  Degradation 
T ransfer  Function 


Zeros  of  the  tranfer  function 
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Incorrect  Knowlegde  of  the  PSF 

POSITIONING  ERRORS  -  25  PT  PSF 


10  pts  w /  incorrect  positions 
Ox  =  Oy  =  4  pixels 


max{an}  =  252 


Effects  of  Impulse  Response  Errors 


Summary  Remarks 


1)  The  zeros  of  the  transfer  function  occur  at 
isolated  points  --  their  effects  on  reconstruction 
are  greatly  reduced. 

2)  Knowledge  of  the  locations  of  the  impulses 
is  critical  to  the  reconstruction. 

3)  Knowledge  of  the  magnitudes  of  the  impulses 
is  not  critical  to  the  recovery. 
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Our  research  is  concerned  with  investigating  the  recovery  of 
intensity  images  that  have  been  degraded  by  a  point  spread  function  (psf) 
consisting  of  a  set  of  irregularly  spaced  impulses.  Our  interest  in  this 
problem  two-fold: 

A)  There  are  some  interesting  aspects  involved  with  characterizing  the 
transfer  function.  In  particular  we  are  concerned  about  the  number 
and  distribution  of  the  zeros  in  the  transfer  function. 

B)  We  feel  that  this  problem  may  have  application  to  the  problem  of 
imaging  through  turbulence. 
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We  are  assuming  a  linear  intensity  based  imaging  system  where  the 
observed  intensity  image,  g(x,y),  is  modeled  as  the  convolution  of  the 
original  scene  intensity  image,  f(x,y),  with  an  intensity  psf,  h(x,y),  plus  an 
additive  noise  term,  n(x,y).  For  our  problem,  we  are  considering  a 
particular  psf  consisting  of  irregularly  spaced  impulses  placed  at  (xn,yn). 
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To  illustrate  this  problem  we  have  chosen  a  picture  of  a  humming 
bird  because  of  the  details  present  in  the  wing  and  neck  feathers  as  well  as 
the  low  frequency  regions  in  the  background.  A  typical  psf  is  shown  here 
with  25  pts  placed  within  a  64  x  64  pixel  region  where  the  image  size  is 
256  x  256  pixels.  The  resulting  blurred  image  is  shown  in  the  lower  left 
while  the  recovered  image  is  in  the  lower  right. 
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Our  primary  interest  in  this  is  to  characterize  the  transfer  function. 
One  particular  area  of  interest  is  to  study  the  number  and  distribution  of 
zeros  in  the  transfer  function.  If  the  transfer  function,  H(fx,fy)  is  rewritten 
in  terms  of  the  real  part,  R(fx,fy),  and  the  imaginary  part,  X(fx,fy),  where  R 
and  X  are  real  functions,  then  in  order  for  H(fx,fy)  to  be  zero  both  R(fx,fy) 
and  X(fx,fy)  must  be  simultaneously  zero. 
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The  real  part  and  the  imaginary  part  are  three-dimensional  functions 
and  therefore  we  expect  the  zeros  to  form  a  set  of  closed-contours.  These 
contours  are  shown  here.  The  zeros  of  the  transfer  function  are  given  by 
the  intersection  of  these  two  sets  of  curves  and  thus,  in  general  are 
expected  to  be  isolated  points. 
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We  are  also  interested  in  the  effects  of  various  errors  in  our 
i  knowledge  of  the  degadation  psf.  One  such  error  is  due  to  incorrect 

positioning  of  the  impulses.  Errors  (An,Pn)  were  applied  to  the  impulses 

where  An  =  fV,  =  0  and  o*  =a^  =4  pixel  . 
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The  effects  of  these  position  errors  can  be  seen  when  the  degradation 
|  psf  is  plotted  along  with  the  psf  used  in  the  recovery. 
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'  Another  source  of  error  can  arise  in  measuring  the  magnitude  of  the 

^  impulses.  To  simulate  such  errors,  the  magnitudes  of  the  recovery 

t  impulses  { an }  were  set  such  that  an=125  and  a3n=75,  resulting  in 

i  magnitudes  in  the  range  10  to  252. 
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The  effects  of  these  two  sources  of  error  can  be  seen  in  the  recovered 
images.  Clearly,  having  correct  knowledge  of  the  positions  of  the  impulses 
1  is  more  critical  than  knowledge  of  the  amplitudes. 

^  Page  10  (Summary) 

1)  The  zeros  in  the  transfer  function  of  the  multiple  point  psf  occur 

i  in  general  at  isolated  points  and  therefore  present  little  difficulty  in  the 

I  reconstruction  problem. 

2)  Knowledge  of  the  correct  positions  of  the  impulses  is  more  critical 
than  knowledge  of  their  magnitudes. 
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